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Abstract 

We study numerically the phase structure of the Ginzburg-Landau model, 
with particular emphasis on mass measurements. There is no local gauge in- 
variant order parameter, but we find that there is a phase transition charac- 
terized by a vanishing photon mass. For type I superconductors the transition 
is of 1st order. For type II 1st order is excluded by susceptibility analysis, 
but the photon correlation length suggests 2nd order critical behaviour with 
z/ ~ ^. The scalar mass, in contrast, does not show clear critical behaviour in 

the type II regime for V ^ 00, contrary to the conventional picture. 
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I. INTRODUCTION 



One of the most interesting phase transitions known is that some materials become super- 
conductive at low temperatures. While the microscopic dynamics behind the phenomenon 
is complicated, there exists a simple effective theory for describing this transition. The 
effective theory is just the Ginzburg-Landau (G-L) model, or a three-dimensional U(l) + 
complex scalar gauge theory. The modulus squared of the scalar field represents the density 
of superconductive electrons. The G-L model might have other applications as well, such as 
the nematic-to-smectic-A transition in liquid crystalsB. 

There exists an extensive literature on the phase diagram of the G-L mo Defining 
the standard G-L parameter as x = m|^/2m^ (where mn is the inverse scalar correlation 
length, or coherence length, and rriw is the inverse vector correlation length, or penetration 
depth), the conventional picture is that the transition is of first order for small x (type I su- 
perconductors), gets weaker with increasing x, has a tricritical point at some Xc, and remains 
of second order for x > Xc (type II superconductors)!. This picture is based on perturbative, 
renormalization group and lattice studies, often in a dual theory. These studies have nev- 
ertheless not been conclusive and there have been arguments concerning, for instance, the 
universality class of the provisional second order transition in the type II regime {x > 1/2)1. 

The purpose of this paper is to study the phase diagram numerically. We improve sig- 
nificantly upon earlier numerical resultsi'i, by having a much finer lattice (smaller lattice 
constant a in physical units) and by measuring, for the first time, the different correlation 
lengths. While the infinite volume and continuum extrapolations are numerically demanding 
and thus the conclusions based on a series of finite lattices can never be quite complete, we 
nevertheless find indications of quite an unexpected pattern in the type II regime. We also 
point out directions for future investigations of this issue. 

The paper is organized as follows. In Sec. || we specify the model studied, in Sec. |Tl| we 
discuss how it is discretized, and in Sec. |I^ how the simulations are organized and what our 
results are. Sec. [V] is a discussion. 
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II. THE MODEL 



Let us first define the model unambiguously. It is a locally gauge invariant 3-dimensional 
continuum U(l) + complex scalar field theory defined by the functional integral 



Z = J VA,V<j)exp[-S{A,<j))], (1) 
S= f d'x\^{diA^ - d,A,)^ +m + tesA,)<j)\^ + 

(2) 



+ml(j)*(f) + A3 (0*0) 

The couplings e|, A3 have the dimension of mass (in units h = c = 1) and, factoring out one 
scale (e|), the free energy density / of the model depends on the two massless ratios 

7/ = x = ^, (3) 

so that 

Z = exp[-Velfiy,x)]. (4) 

The phase diagram of the theory can thus be drawn in the {x, y)-plane. 

Since the theory in eq. (Q) is a continuum field theory, one has to consider ultraviolet 
renormalization. There is a linear 1-loop and logarithmic 2-loop divergenc for the mass 
parameter m|. In the MS dimensional regularization scheme in 3 — 2e dimensions, the 
renormalized mass parameter becomes^ 

2 -4e| + 8A3e^ - SAj ^ 
mM = log—, (5) 

where specifies the theory. To be more precise, we thus define y = ml{el)/e^, which 
specifies the continuum theory at the full quantum level equally well as A^,. 

The physical values of x and y depend on the microscopic theory behind the effective 
theory in eq. (j^). For reference, for usual superconductors in the notation of Ref.i, 

y=—(L_l] ^ = ^^^ (6) 
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FIG. 1. The qualitative phase diagram of the G-L theory. 

with g ~ 10^^, g ~ 0.01, r < 1. For high superconductors x can be ^ 1. At present we 
consider the theory (g) as such and questions of vahdity, hke the need to include (0*0)^- 
terms, do not enter. 

The phase diagram of the G-L model on the y, x > plane (see Fig.|l]) contains a curve 
y = yc{x) along which the system has a first order transition for small x. To one vector 
loop this curve is simply given by yc = l/(187r^x). We shall call the region y > ydx) the 
normal and the region y < ydx) the superconducting (SC) phase. At small x the reliability 
of perturbation theory has been verified numericallyEl. However, perturbation theory gets 
worse at large x and the issue now is what happens then. 



III. DISCRETIZATION 



To latticize (a = lattice spacing) the theory with the fixed continuum variables x, y, we 
introduce the link field f/j(x) = exp[me3Aj(x)] = exp[iaj(x)]. Relating the counterterms in 
the MS and lattice regularization schemeJi^, the lattice action be 



jecomesi^ 



S = (3gT. [1-cos4(x) 

X,i<j 



-/5j,5lRe0*(x)[/,(x)0(x + 



X,i 



y 3.1759115(1 + 2a;) 
Wg 

-4 + 8a; - 8a;2)(log6/3G + 0.09) + 25.5 + 4.6a;- 



^ X 



167r2/32 



+ ^E[r(x)0(x)]^ (7) 
4Pg X 

where /3h is so far arbitrary, /?g = l/e|a, and 

Fij(x) = aj(x) + aj(x + i) - ai(x + j) - aj(x). (8) 

We have scaled the continuum scalar field to a dimensionless lattice field by (j)*(j) — > 
PH(p*<P/'^(^y but further rescalings (by a specific choice of Ph) are possible; we thus scale 
the coefficient of 0*0 to be +1. For a given continuum theory (e|,?/, x), eqs. (0)-(|D spec- 
ify up to terms of order e|a the corresponding lattice action. It should be noted that the 
complicated counterterm expression in the square brackets in eq. (^ only affects the value 
of Uc for given x (the counterterm guarantees that the limit yc{a — >• 0) exists), but not the 
qualitative structure of the phase diagram. 

One of the most essential points of the present lattice simulations is the extrapolation 
to the continuum limit: first the infinite volume ^ oo at fixed lattice spacing a, then 
a — > 0. To estimate the required sizes of V = (Na)^ and a, consider a system with a typical 
correlation length ^. Then one has to satisfy (on a periodic lattice) a <^ ^ Na/2 or, in 
physical units, 

e^a = — < el^ < -— . (9) 

We observe that ejC, ~ 1 and take /3g = 4, 6 so that 1//3g <^ ej^; and = 32, . . . , 64 so 
that N/{2(3g) > eli. Note that RefB had (3g = 0.2, < 15, so that the lattice spacing a 
was larger than the typical correlation lengths. 

Apart from the UV-cutoff effects discussed, there is another effect related to a finite a. 
Indeed, we use a compact formulation for the U(l) gauge field, which changes the topol- 



ogy of the continuum theory and imphes that the photon becomes massive0. However, 
semiclassical computation for this Polyakov mass giveslll 



63 



exp 



3.1767r 



4 



G 



(10) 



Thus for /3g ^ 4 this photon mass (as well as the monopole density from which it origi- 
nates) should be negligible {m^ /e\ < 0.01) and our results the same as in the non-compact 
formulation within statistical errors. 



IV. SIMULATIONS AND RESULTS 

For the simulations we choose two values oi x = 0.0463 and x = 2, corresponding to 
a strongly type I and type II superconductor. We then measure averages of local or bilocal 
gauge invariant quantities and locate the critical curve y = ydx) on which the system 
changes its properties. Note that there is no local gauge invariant order parameter which 
would vanish in either of the phases. However, we shall find that the photon mass, measured 
from a correlator, vanishes in one of the phases within errors. 

The phase transition is located by finding the maximum in y of the susceptibility x defined 

by 

X = elv((^-(^)Y), (11) 

where 0*0 is the volume average, 0*0 = Jd^x(j)*(j), and by studying its large- 1/ behav- 
10 r@. If there is a first order transition, the distribution of 0*0 precisely at ydx) has two 
peaks which remain at fixed distance and get narrower when V ^ 00. Then the maximum 
of X grows as the volume V. In a second order transition the expected behaviour is ~ V^, 
K being a critical exponent. If x ~ then either k < 0, or the transition is of higher than 
second order or absent. 

The susceptibility maximum is plotted in Fig. |^. One sees a very clear difference between 
X = 0.0463 and x = 2. The behavior of the system at x = 0.0463 at the largest volumes 
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FIG. 2. The maximum of the susceptibihty x as a function of volume. The straight Une is ~ V. 

indicates a linear 1st order behavior. At x = 2, in contrast, the transition is not of first 
order. If the transition is of second order the critical exponent k, is close to zero, as noticed 
already ini'§. However, a still higher order transition (or a smooth crossover) cannot be 
excluded based on these measurements. We thus turn to masses. 

For the mass measurements we use the scalar operator (f)*(f) and the two vector operators 
(f)*Di(j) and eij^Fj^. On the lattice these are 

O(x)=0*(x)0(x), (12) 
0,(x) =Im0*(x)f/,(x)0(x + 2), (13) 
di{x) = eijk sin Fjk{x). (14) 

The correlation masses are measured (from a lattice of size N^N^) by first summing over 
planes, possibly with momentum p = 27r/(aA^^), 

0{z;p) = Y.e'^^0{x,y,z), (15) 

x,y 

and then studying the large-t behavior of 



G{t; p) = -i- Y.{0{z; p)0*{z + t;p)). (16) 



The momentum is needed for the correlator of O3, used to measure the photon masgi^l (it can 
also be used for the very light scalar mass, to get a better signal in a finite volume): without 
the factor e*^^ in eq. (|15|), the plane average for O3 would simply vanish. In perturbation 
theory, one finds for the asymptotic behavior of this correlator in the continuum limit. 



where = in the normal phase. At 1-loop level. 



A^ = l ^- — . (18) 

Very close to the critical point ?/ ~ the expansion for thus breaks down, and the form 
of the correlation function is determined by an anomalous dimension. For our datapoints, 
the functional form in eq. (|T^ fits the data well and the energy E is measured from the 
exponential fall-off. 

To improve the projection to the low-lying mass states, it is indispensable^ to use blocking 
techniquesii to define extended operators, and to make a mixing analysisH between operators 
at different blocking levels, to search for the linear combination giving the best signal. We 
find the best results with blocking level 3 for the scalar mass. The mixing analysis allows to 
get a good signal already at a relatively small t, t ~ (2 — 5)a. We denote by m-y the lowest 
lying vector state, and by mw the first excitation (where it can be determined). 

The masses (= 1/^) near the transition are shown in Fig. ^. For type I superconductors 
(x = 0.0463) one observes in the SC phase one scalar with a rather small mass mu- The 
vector operators couple dominantly to a single state of a larger mass mw = ttl^. These are 
the standard (inverse) coherence length and penetration depth. Close to the transition point 
one can observe both phases simultaneously, and the masses are discontinuous. Above the 
transition the scalar excitation is there with a larger mass, but the vector operators couple 
to two quite different states. There is a 0*0 bound state of large mass rriyy, seen dominantly 
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by the operator Oj, while the operator O3 sees the photon of mass which is consistent 
with zero within 2a. The picture here is the standard one of a first order transition with 
as an effective order parameter. 

For type II superconductors (a; = 2) a 1st order transition was excluded by the suscep- 
tibility analysis: no two-peak structure exists. However, Fig. |^ shows that there is some 
transition since in the normal phase = within errors. The critical region is shown 
in a magnified form in Fig. |, and one can see that could go to zero continuously (on 
the other hand, it should be noted that a discontinuity in cannot be excluded either). 
The dashed curve in Fig. ^ shows the fit m^/e\ = A{yc — y)" over the whole y-range, where 
A = 0.95(2), yc = —0.046(6), u = 0.44(2). For the lower solid curve u has been fixed to 
z/ = In Fig. ^ for a smaller y-range, v = 0.39(17). The behaviour of is thus consistent 
with a mean field exponent. 

As for rriw, it is seen that it deviates from already below the critical point, unlike in 
the first order case shown in the top panel of Fig. |^. This might indicate that the transition 
line has split into several transitions in the type II regime. However, mw is an excited state 
and thus there are some systematic errors in its determination which may be larger than 
the statistical errors shown in the figure. 

Consider then the scalar mass rriH- First, note that its minimum is at a point different 
from where diverges (see Fig. Second, rriH dips steeply in Fig. |], but does not go to 
zero as can be seen in Fig. |[ On a finite (periodic) lattice with spatial extension L = Na 
one cannot expect to see scalar mass values smaller than ~ 2/L. However, according to 
Fig. ^ the scalar mass is larger than this value and does not show the corresponding volume 
dependence. In fact, mn is volume-independent within error bars at the minimum. 

The question remaining then is whether the scalar mass depends on the lattice spacing a. 
This seems unlikely, since the lattice spacing we used is much smaller than the correlation 
length, e|a = 1/Pg = 1/4 << e|^jy = e^/rriH ^ 2, so that one does not expect large effects 
from removing the UV cutoff. Indeed, we have made simulations with /?g = 3 and /3g = 6 
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at a few points around the minimum, and we do not find any appreciable lattice spacing 
dependence. It should be pointed out, though, that for /^g = 6 one should go to ~50% 
larger lattice sizes N to get physical volumes comparable with those for jSc = 4, according 
toeq. i). 

Note also that as long as a 7^ there is in a strict sense no phase transition due to the 
Polyakov mass. However, the Polyakov mass in eq. (|To|) is very small for = 4, so that its 
finite value should have no effect. 

V. DISCUSSION 

We have seen that in the type II regime, the data are consistent with a second order phase 
transition driven by a diverging photon correlation length. In contrast, the scalar mass 
shows quite unexpected behaviour. It appears that rriH has a minimum away from where 
the photon correlation length diverges, and that it does not show critical scaling towards 
zero when the volume increases. If this pattern remains there for larger volumes and smaller 
lattice spacings, then one has to modify the standard picture of the superconducting phase 
transition in the type II regime. 

As to the critical exponents of the transition, we have measured the finite size scaling 
susceptibility exponent defined after eq. ([TT| ) to be consistent with zero, and the photon 
correlation length exponent in the SC phase to be consistent with |. We have not measured 
the anomalous dimension 7]a of the photon correlation length, (O3(x)O3(0)) ~ l/|x|^+''^, 
but in principle it could be measured with a significantly extended analysis (it would show 
up in the functional form of eq. (0)). If the transition were of 2nd order in the usual sense, 
then one could also measure the scalar correlation length critical exponents u, v' in the two 
phases, mn ~ |y — Z/cT''^', and the corresponding anomalous dimension at yc- However, as 
we have not seen any critical behaviour for the scalar correlation length, these exponents 
cannot be systematically measured. It seems that away from the critical point, mn shows 
approximate scahng where i/, z/' are consistent with |, see Fig. |^. 
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A more conclusive study of the infinite volume and continuum limits in the type II regime 
would clearly be needed. Unfortunately this limit is numerically very demanding and re- 
quires much more extensive further simulations. Based on the present investigation, we can 
nevertheless point out that it would probably be more economic to use the non-compact 
lattice action than the compact one we have used here. The reason is that one can then 
take a large lattice spacing (smaller Pg) in order to get a larger physical volume, without 
having to worry about the non-zero value of the Polyakov mass. 

As a final observation, let us note that one might expect the study of a gauge theory 
with a U(l) gauge group to be simpler than that with a more complicated group such as 
SU(2)i0 or SU(2)xU(l)ii. In fact, this turns out not to be the case: U(l) is numerically 
more demanding and requires larger lattices (at least in the compact formulation). The 
existence of a massless photon is not as such the only reason: there is a photon also in 
SU(2)xU(l), but there it exists in both phases and is not an order parameter. 
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FIG. 3. The masses at x = 0.0463 (top) and x = 2 (bottom). Note the discontinuous and 
continuous behaviour of my in the top and bottom figures, respectively. The values of m-y in 
the top figure for y = 0.12, 0.14 are 2. . . 3 standard deviations from zero, but all or most of this 
deviation is expected to be due to the absence of a mixing analysis at these data points. The curves 
in the bottom figure represent fits ~ A\y — yd'^. The dashed curve for rrij is with a free exponent 
u = 0.44(2) over the whole y-range whereas for the solid curves, u has been fixed to u = ^. 
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FIG. 4. A magnification of tlie region around tlie critical point for x = 2 in Fig. |3[ The dashed 
curve is a fit with a free exponent over the range y = —0.25. ..0.0, and gives = 0.39(17). The 
sohd curves are with v = ^, and with a constant for itih (see Fig. 
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FIG. 5. The volume dependence of the scalar mass. The results are consistent within er- 
ror bars at the minimum. The expected critical scaling behaviour is shown with the horizon- 
tal lines indicating the value of 2/L. No such scaling is observed. The fits are of the form 
niH/el - U'(yc - y)^/^ + b] {y < y^), U(y - yc)^/^ + b] {y > y^), with A' - 1.1, A - 1.7. 
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